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Abstract 

The propagation of an intense femtosecond pulse in a Raman-active medium is analyzed. 
An analytic solution which describes in explisit form the evolution of the light pulse is 
derived. The field of an intense light wave undergoes a substantial transformation as the 
wave propagates through the medium. The nature of this transformation can change over 
time scales comparable to the period of the optical oscillations. As a result, the pulse 
of sufficiently high energy divides into stretched and compressed domains where the field 
decreases and increases respectively. 


1 Introduction 

The physics of interaction of intense ultrashort light pulses with nonlinear media has attracted 
interest because of progress in subpicosecond-range laser technology and the attainment of laser- 
beam power levels of terrawatt range (see, for example, [1]). The light pulse dynamics and the 
medium evolution in an intense field differ qualitatively from the usual picture drawn by standard 
nonlinear-optics perturbation theory. A fundamental distinctive feature of ultrashort pulses is 
that their duration is shorter than the time scale of the response of the medium, so the interaction 
definitely occurs in a coherent regime. The standard approximation of a slowly varying amplitude 
and a slowly varying phase of the field becomes ineffective. A description of the interaction based 
on the actual (instantaneous) field values is appropriate [2]. 

The Raman-active media can be effectively excited by a single femtosecond pulse because its 
wide spectrum initially contains intense Stokes and anti-Stokes components of the field [3, 2]. The 
qualitatively new interaction regime, so-called self-scattering, occurs for light intensities much 
lower than the threshold ones for ordinary stimulated Raman scattering developing from sponta- 
neous noise [4]. The estimations testify to the fact that the effective excitation of high vibrational 
levels and even the dissociation of molecules in the field of an ultrashort electromagnetic pulse 
can be observed using modern femtosecond lasers [5]. 

On the other hand, the shape and spectrum of femtosecond pulse should undergo a substan- 
tial transformation as it propagates through the Raman-active medium. The description of pulse 
evolution should take into account a substantial redistribution of the medium level populations 
during the pulse duration that results in different interaction regimes for different pulse fractions: 
either absorption due to Stokes scattering or amplification due to anti-Stokes component genera- 
tion become predominant. The simplest model which is widely used in the theory of stimulated 
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Raman scattering and includes the effect of medium saturation is the model of two-level nonlin- 
ear oscillator. Within the framework of this model we succeeded to describe in explicit form the 
evolution of high-power femtosecond pulse in Raman-active medium. 


2 Solution of Wave Equation 


We shall describe the dynamics of a Raman-active medium in the field of an ultrashort pulse by 
a two-level model of a nonlinear oscillator [6]: 


d 2 Q 

dt 2 


+ _L = — L 

T 2 dt W 2 M 




dp P~ Po 
dt T\ 



(1) 


Here Q is a normal coordinate, M is the reduced mass, fl is the eigenfrequency of the equivalent 
nonlinear oscillator (the Stokes shift), Ti and T 2 are relaxation times, the coefficient da/dQ is the 
derivative of the polarizability at the equilibrium value Q = Q 0 , and p is the difference between 
the populations of the upper and lower levels (the value p 0 = p(t — — oo) corresponds to the state 
of the system before the beginning of the interaction with the field). 

Equations (1) are to be solved jointly with the wave equation. In the case at hand, in which 
the interaction of the carrier pulse with only the scattered wave propagating in the same direction 
is taken into account, the wave equation can be written 

d£ \_dE_ _ _ 4tt dP_ 
dz + c dt c dt 

with macroscopic nonlinear polarization of the medium P = N EQ induced by the field 
E(z,t), where N is the density of the medium. We stress that the quantity E — E(z,t) in (1) 
is the instantaneous value of the pulse field strength, not its envelope. The reduced form of the 
wave equation (2) is provided by the fact that the stimulated Raman backward scattering is very 
weak owing to the short interaction length of pulses moving in different directions [6]. 

We restrict the analysis to the case of the coherent interaction, in which the pulse duration 
does not exceed the time scales of the response and relaxation of the medium r p <C ft -1 , Ti, T 2 . 1 

Material equations (1) can then be integrated for an arbitrary time dependence of the field 
E(z , t ) and the solution looks like rotation of material variables 

... dQ ( hCl Y' /2 . /da\ ... . 

p(z,t) = ft, cos ♦(z, I), = -po ( sign sin (3) 


The phase ty(z,t) of the material variable rotation is directly proportional to the energy of the 
pulse fraction which has passed through the given space point z up to the given time t 

(2 hQM)- 1/2 f E 2 {z,t , )dt ' (4) 

'Strictly speaking the coherent regime of interaction is provided by the condition Qo = | |(2fif2M) -1 / 2 £'o > 

Q . , T, - 1 , T 2 _ 1 , i.e. , the analog of Rabi frequency of the two-level oscillator is the highest frequency of the problem. 
For high-energy femtosecond pulses this condition as well as inequality r p <C f2 -1 , T\, T 2 leads to the solution 

(3), (4). 
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dQ 
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From (3) it follows a strong nonlinear dependence of polarization on the pulse field and the wave 
equation of the form 


dE 

dz 


1 dE 

c dt 


1 + fie f sin ty(z,t , )dt > 

J — oo 


— (3E sin 'k, 


(5) 


where | N | da/dQ \ ( hQ/2Mc 2 ) 1 / 2 is the inverse length of induced Raman self-scattering. 

The nonlinear equation (5) allows an analytic solution describing in explicit form the evolution 
of the pulse shape and spectrum [7]. 

Multiplying (5) by E and integrating it with respect to time one can easely obtain the equation 
for phase ty(z,t) (i.e. , for current pulse energy) 


1 d ^ 

dz c dt 


1 -f 0c f sin t f ) dt ' 

J — oo 


= /?(cos — 1 ). 


( 6 ) 


Let us regard the pulse field E(z,t) as a function of the spacial variable 2 and the phase 'I'( 2 ,<), 
i.e., E(z,t) = E(z,ty). Taking into at 
rewrite Eq. (5) in the following form: 


i.e., E(z,t) = E(z , $). Taking into account that |f = §f + §f § 7 , = ff fp, we can finaly 


3E a/ ,, 8E 

— + / J(cos*-l)^ = 


-f3E sin 'P. 


(7) 


This partial differential equation can be easily solved by integrating along characteristics on 
which 
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where is the given phase at the boundary of the medium (i.e., the energy of the pulse fraction 

which has entered the medium) which is connected with the field strength at the boundary by the 


relation ^o(^) = 
corresponds to the p 


( 2 /jfIM ) -1 / 2 Eq(t]') dr}'. According to ( 8 ) the pulse fraction which energy 

ihase 'to = 2?rn, n = 1,2,... moves through the Raman-active medium under 
conditions of self-indused transparency, when the energy of the leading part of 27r-pulse absorbed 
due to Stokes scattering completely returns to the trailing part of the pulse due to anti-Stokes 
scattering. In this case the spectrum of the leading part of the pulse becomes enriched by long- 
wavelength components of the field and the spectrum of the trailing part - by the short-wavelength 
components. This analysis generalizes the results found on the basis of numerical calculations [ 2 ], 
The characteristics themselves are given implicitly by the expressions 


7 / -|- $( 2 , 7 /) = t — zjc , ( k( 2 .ri'| = dz 


f drf 
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Q ^ 
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( 10 ) 


where $(2,7/) is the nonlinear delay of the individual parts of the pulse as it moves a-way from the 
boundary. 
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It is convenient to describe the evolution of the pulse shape by an effective frequency u>{z,t) 
of the pulse field oscillations. The value u>(z , t ) characterizes the density of field oscillations for 
different parts of the pulse and its variation as the pulse propagates through the medium. It 
follows from (10) that if we fix a small fraction of the pulse its duration At at the space point 2 is 
connected with that one Ar 0 at the boundary by the relation At = A t 0 (1 + d$/dr /) . The value 
At determines the field oscillation period at the given space point 2 , hence the effective frequency 
of the field oscillations transforms in accordance with 


u>(z,t) = 


1 +0Z 


Mv) 

Bz 

sin ^0 + — (1 - cos ^ 0 ) 


( 11 ) 


3 Femtosecond Pulse Evolution in Raman-active Medium 

It can be seen from Eqs. (9) and (11) that the changes in the field E(z,t) and the characteristic 
frequency u>(z,t) over space occur identically. At the beginning of the pulse, when the phase 
of the two-level oscillator satisfies « 1 (be., when the energy of the pulse fraction which 
has entered the medium is small), there are decreases in the field strength and the oscillation 
frequency: E(z,t) = Eq(tj) ( 1 + /3z^ 0 /2) -2 , u(z,t) = Uo{v) (1 + f3z ) $o/2)~ 2 . This case corre- 
sponds to Stocks scaterring. Later, when the phase becomes greater than 7r and reaches the value 
arctan(— 2 / fiz), the field and its frequency increase. The generation of anti-Stokes components of 
the field thus becomes predominant. At the nature of the pulse transformation changes 

again. At a given point in space at different times, corresponding to different characteristics (10), 
we thus observe oscillations of regimes of compression and stretching of the field oscillation peri- 
ods. The sequence of regimes of stretching and compression of the pulse with increasing value of 
the incident energy reverses when we switch from an originally absorbing medium to an originally 
inverted one. 



FIG. 1. 
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A descriptive way to analyze the pulse evolution is provided by Fig. 1. It is seen from (9) 
and (11) that the nature of pulse transformation is determined by the sign of the expression in 

brackets of the denominator sin + ^r(l — cos 'ko) = 2 sin 2 ^ |(tan . If 'ko < tt, we 

have (tan > — ^ for any distance 2 from the medium boundary. It means that the energy 

of such pulses can be only absorbed and the absorption is accompanied by the increase of the field 
oscillation periods for any pulse fraction (and as a result, by the increase of pulse duration) and 
by the shift of the pulse spectrum into the low-frequency region due to Stokes scattering. 

If 7 r < $0 < 27 r, there are two different regimes of pulse transformation. The leading pulse 
fraction 0 < 'ko < < 4 /q = arctan (— 2/^z) undergoes energy absorption and stretching of the 
field oscillation periods and its spectrum becomes enriched with low-frequency components of the 
field. The field in the trailing pulse fraction amplifies during the pulse propagation, the oscillation 
periods decrease (that results in the compression of this pulse fraction) and the spectrum is 
enriched with high-frequency components. For 27r-pulse all the energy concentrates in the trailing 
edge which is compressed as the pulse propagates through the medium. The same picture takes 
place for pulse fractions 27rn < 'ko < 2n(n + 1). 

Figure 2 shows an example of pulse evolution described by the solution (8)-( 1 1 ) for 'ko = 47r. 
The input pulse represents a two-period fraction of a sine-shaped signal. Figure illustrates the 
time dependence of the pulse field at different distances from the medium boundary. Thus, the 
high-energy pulse is devided into the set of compressed powerfull 27r-subpulses the maximum field 
strength and inverse width of which at large distances increase ~ (fdz) 2 according to (8)-(10). 
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